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Charge transport in bacteriorhodopsin monolayers: The contribu- 
tion of conformational change to current-voltage characteristics 



E. Alfinito^'^ and L. Reggiani^'^ 



ON 
O 
O 
(N 

u 
en 



d 

•1— I 

(N 
> 

O 
(N 
l> 

o 

oo 

O 



% 



^ Dipartimento di Ingegneria delVInnovazione, Universita del Salento - via Arnesano, 1-73100 Lecce, Italy, EU 

^ CNISM - Consorzio Nazionale Interuniversitario per le Scienze Fisiche della Materia - via della Vasca Navale 



i,I-00146 Roma, Italy, EU. 



PACS 87 . 85 . j c - Electrical, thermal, and mechanical properties of biological matter 
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Abstract. - When moving from native to light activated bacteriorhodopsin, modification of 
charge transport consisting of an increase of conductance is correlated to the protein confor- 
mational change. A theoretical model based on a map of the protein tertiary structure into a 
resistor network is implemented to account for a sequential tunneling mechanism of charge trans- 
fer through neighbouring amino acids. The model is validated by comparison with current-voltage 
experiments. The predictability of the model is further tested on bovine rhodopsin, a G-protein 
coupled receptor (GPCR) also sensitive to light. In this case, results show an opposite behaviour 
with a decrease of conductance in the presence of light. 



Introduction. — In recent years, several papers re- 
ported on the electrical properties of single proteins in- 
serted into hybrid systems [1-5]. The principal aim of 
these researches is to explore the possibility to design nan- 
odevices, mainly nanobiosensors, with extreme sensibil- 
ity and specificity. Most work was devoted to metallo- 
proteins, which exhibit semiconductor-like conductivities, 
but increasing attention has been also addressed to bac- 
teriorhodopsin (bR), the light sensitive protein present in 
Archaea. Bacteriorhodopsin exhibits a conductivity close 
to that of an insulator, and is robust against thermal, 
chemical and photochemical degradation [3,4]. It is a pro- 
tein (opsin)-chromophore complex that underlies photon- 
activated modifications in which the chromophore changes 
its structure from the all-trans to the 13-cis. The chro- 
mophore evolution induces multiple protein transforma- 
tions, which fastly drive the protein from the Ksgg to the 
L550 state. In the slow transition from this state to the 
M410, a proton is released. Junctions prepared with bR 
monolayers contacted with Au electrodes, show non-linear 
/- V characteristics both in the dark (or blue light) and in 
the presence of green light [3-5] . It is generally assumed 
that green light produces a more intense current as a 
consequence of the modification of the chromophore-opsin 
complex. Although the experimental results are not com- 
pletely assessed from the quantitative side [4], the main 



results of available experiments can be summarized as fol- 
lows, (i) A conductivity about four orders of magnitude 
higher than that of a homogeneous layer of similar dielec- 
tric material was observed, (ii) The current voltage (/- V) 
characteristics exhibit a super-Ohmic behaviour pointing 
to a tunnelling charge-transfer mechanism, (iii) The pres- 
ence of green light significantly increases the electrical con- 
ductivity up to about a factor of three with respect to the 
dark value, (iv) The I-V sensitivity to the green light 
is washed out by substituting the original chromophore 
with one that is not sensible to light, (v) When the chro- 
mophore of the protein is completely removed, the conduc- 
tivity is suppressed by three orders of magnitude lower, 
taking the value of a standard dielectric. Despite the 
presence of these experiments, the existing theoretical ap- 
proaches only provide a phenomenological interpretation 
of some I-V characteristics in terms of metal-insulator- 
metal tunneling theories [5,6]. 

The aim of this work is to develop a theoretical model 
able to capture the main features of the above experiments 
and sufficiently general to constitute a new framework 
for describing electron transport properties through pro- 
teins. To this purpose the protein conductivity is modelled 
by a sequential tunneling mechanism through the protein 
amino acids, whose positions pertain to a specific confor- 
mation. Thus, the conductivity is directly connected to 
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the protein tertiary structure and a light induced struc- 
tural change will result in conductivity change which can 
be quantitatively determined. 

The network tunneling model. As general fea- 
tures, the protein electrical conductivity is close to that 
of an insulator, but it can increase significantly in par- 
ticular environmental conditions, like in organic semicon- 
ductors [7]. Furthermore, the microscopic mechanism of 
charge transport in proteins is very articulated and in- 
cludes electron, ion, and proton transport [8] in a yet 
not completely understood sequence. The difficulty to 
discriminate among these kinds of transport is further 
enhanced by the lack of experiments, mainly related to 
the difficulty of working with biological materials at the 
nanoscale. The most qualified mechanisms to describe 
electron transport in proteins are the hopping/tunneling 
processes [9,10]. In this framework, the transport model 
which assumes the existence of privileged pathways of se- 
quential tunneling has gained a wide consensus [11]. By 
embracing this picture, we describe the protein as a net- 
work and the electron transport as due to a sequential tun- 
neling along preferential channels. In doing so, we produce 
a map of the protein by taking from the protein data base 
(PDB) [12] its tertiary structure, and setting up a graph 
whose nodes correspond to the amino acids (identified by 
the a-carbons) [13,14]. Two nodes are then connected 
with a link only if they can interact within a given inter- 
action radius Re- 

At this stage, we neglect thermal fluctuations because 
the relevant signals of experiments do not suggest their 
evidence. 

To the purpose of investigating the conduction prop- 
erties of the protein, the graph is turned into a resis- 
tor network: Each link is replaced by an elemental resis- 
tance, which mimics the different charge transfer property 
between neighbouring amino acids associated with their 
mutual distance. To detect the differences in conduction 
properties between two protein states, as those produced 
by the conformational change, we adopt the simplest pa- 
rameter model by choosing the same resistivity for each 
amino acids. In this way, the tertiary structure and the 
interaction radius arc the relevant input parameters of the 
theory. 

As elemental resistance of the network we take 
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where Aij ~ '^{R'i ^ ^j/'^) is the cross-sectional area 
between the spheres of radius Re centered on the i-th and 
j'-the node, respectively, Uj is the distance between these 
centers, p is the resistivity. 

By construction, each elemental resistance depends 
upon the distance between nodes. Therefore, following up 
a conformational change, the variation of this distance im- 
plies a variation of each elemental resistance, which even- 
tually leads to a variation of the network resistance (and 



thus of the protein resistance) . As a consequence, a topo- 
logical transformation can be monitored by means of re- 
sistance measurements. 

By taking the first and last amino acid of the protein 
primary structure as ideal electrical contacts, the protein 
resistance for a given Re is calculated by solving the cor- 
responding linear resistance network within a standard 
Kirchhoff framework [14]. Being interested to a change 
of the resistance rather than its absolute value, this choice 
for the contacts does not modify the results of the model. 
Numerical calculations show that for Re below a thresh- 
old value (typically around 4 A ) the network is discon- 
nected and thus not conducting at all. By contrast, for 
Re well above the average distance between amino acids 
(typically above 15 A) [14] each amino acid is connected 
with most of the others, thus the conduction becomes in- 
sensitive to modification of the structure. Between these 
critical values of Re, the change of resistance exhibits a 
smooth behaviour. Accordingly, we select the value of Re 
which maximizes the change of resistance due to the con- 
formational change, and which is found to be around 6 A , 
a well accepted value for the interacting radius between 
amino acids [13-15]. 

To account for the strong super-linear /- V characteris- 
tic, the model implements a barrier-limited current mech- 
anism as follows. At increasing voltage, each elemental 
resistance is allowed to take a second value of resistiv- 
ity which, playing the role of a small series resistance of 
the network, is several order of magnitude lower than the 
first one. The probability of this choice mimics a barrier- 
limited mechanism in analogy with the case of an organic 
molecular layer [6]. In this way, the initial linear increase 
of current with applied voltage turns into a supcrlinear in- 
crease, with the value of the barrier energy to be fitted by 
comparison with experiments. Accordingly, the stochas- 
tic selection is taken to be ruled by the direct tunneling 
probability ^ [6]: 
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where Vi^j is the potential drop between the i-th and j-th 
node, m is an effective electron mass, here taken as that 
of the free electron, and $ is the barrier height. 

We notice that the model does not take into account the 
protein environment mainly for two reasons. The former 
accounts for the fact that here we are interested only in 
the variation of protein response, and this should be inde- 
pendent on the environment, since it remains the same in 
the conformational change. The latter reason is in agree- 
ment with ref. [4], which concludes that transmembrane 
electron transport occurs essentially only via bR and not 
by the lipid bilayer. Furthermore, the chromophore is not 



^We remark that the modeUing leaves the possibihty to replace 
tunneling by a thermionic-emission mechanism, in case further ex- 
periments would evidence a significant temperature dependence of 
transport characteristics. 
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included in the network being its contribution to the con- 
ductivity implicitly considered through the induced defor- 
mation of the structure. Otherwise, we do not exclude the 
possibility that the chromophore can further contribute to 
the total conductivity. 

To make a reliable comparison between the native and 
activated configuration, the corresponding representations 
have to be taken from identical experimental settings [14]. 
Table 1 reports some characteristics of four couples of 
bR representations which satisfy this requirement. Our 
model is strongly dependent on the protein geometry, 
therefore, we focus our attention on the representations 
obtained with the best resolution, (IMOK-IMOM) and 
(2NTU-2NTW). In the former couple IMOK represent the 
K-state instead of the native state, and IMOM refers to the 
M-state. In the latter couple, 2NTU refers to the native 
state, while 2NTW refers to the L-state. Furthermore, we 
have verified that all the four couples reported in table 

1 show similar impedance change, i.e. the native state 
has a larger zero-frequency impedance with respect to the 
activated state, any activation step is considered. In con- 
clusion, we select the couple (2NTU-2NTW) as computa- 
tional sample because it has the best resolution, contains 
the native state and although the activated state is that 
which forthwith comes before the M-state, we know that it 
can only underestimate our results, while preserving their 
substantial validity. 

Calculations proceed as follows. The native and acti- 
vated proteins are mapped into the corresponding topo- 
logical networks. The resistance of the corresponding net- 
work at the given voltage is then calculated using the value 
i?c = 6 A , which is found to give the best contrast. The 
value of p for each elemental resistance is stochastically 
determined on the basis of the probability given in eq. 
(2) between the values pmax and Pmin- The value of 
Pmax = 10"'^'' flm is chosen to reproduce a current of 
about 10^^*^ A per bR monomer at 0.5 V, in accordance 
with experiments [3] . A rescaling of this value simply im- 
plies a rescaling of the total current [6]. Thus, since our 
interest focuses on the variation of current and not on its 
absolute value, the pmax value serves as a normalizing 
parameter. The value of Pmin — 10^ i^m is chosen to 
account for the superlinear behaviour of the /- V charac- 
teristic and serves as series resistance of the circuit. In 
the voltage range considered in the experiments (0-^1 V) 
the deviation from linearity is sufficiently weak so that the 
value chosen for pmin does not play a significant role; in- 
deed, the results do not change even varying its value over 

2 orders of magnitude. On the other hand, we found cru- 
cial the role of the barrier height, which drives the shape 
of the superlinear behaviour. 

The network resistance at the given voltage is then cal- 
culated by using the following iterative procedure. First, 
the network is electrically solved by using the value pmax 
for all the elemental resistances. Second, each pmax is 
stochastically replaced by Pmin using the probability in cq. 
(2) according to the local potential drops calculated in the 



first step; the network is then electrically updated with the 
new distributed values of pMAX.min- Third, the electrical 
update is iterated (typically 10** -^ 10^ iterations depend- 
ing on the value of the applied voltage) by repeating the 
second step until the resistance of the network, taken as 
the average value over the iteration steps, < R >, con- 
verges within an uncertainty less than 1 %. In this third 
step, the initial iterations (100 -^ 2000 depending on volt- 
age) contain a significant numerical noise and as such they 
are disregarded to avoid an unwanted drift of the average 
value (see fig. 1). Finally, the current at the given voltage 
is calculated as / = T^/< i? >. 

The comparison between theory and experiments pro- 
ceed as follows. On the one hand, the theoretical model 
describes the /- V modifications of a single protein when it 
undergoes a conformational change. On the other hand, 
the experimental results are carried out on a macroscopic 
sample of 5 nm width and 2 x 10~^ crn^ cross-sectional 
area. Therefore, to compare the results, we normalize the 
current value of the single protein in its native state at 1 
V to that of the experiment in dark. (The normalization 
corresponds to multiply the current of the single protein 
by a factor 10^ — 10^). The same normalization factor 
is used for the case of the single protein in its activated 
state, which is then compared with the experimental value 
under green light. 

Results. — Figure 1 shows a typical evolution of the 
average resistance of the single protein at low and high 
voltages. Within the required uncertainty, at both volt- 
ages the resistance is found to converge above about 4x 10'^ 
iteration steps, with an initial behaviour noisier at higher 
voltages. This evolution is typically observed for any value 
of $ considered in calculations. For the selected PDB 
entries reported in table 1, calculations give a resistance 
value of the native state larger than that of the activated 
state by about 10 % in the linear regime. 

Figure 2 reports the single protein resistance vs. the 
applied voltage for the native and activated state at two 
different values of the barrier height. Here, we found that 
the sharp drop of the resistance starts at about 0.7 V, 
being more effective for $ = 0.53 meV. Furthermore, in 
the region of voltages of interest for experiments, — 1 V, 
a further decrease of the series resistance is found to play 
a negligible role. 

Figures 3 report a quantitative comparison between the- 
ory and experiments. The experimental results, as de- 
tailed in ref. [4], give the net message that green light 
illumination induces an increasing of current in bR junc- 
tions. The value of this increase depends on the junc- 
tion preparation (bR/PC vescicles or bR/OTG vescicles), 
and has a significant range of variability, from a mini- 
mum of about 30% to a maximum around 300%. Accord- 
ingly, taking into account this experimental uncertainty, 
the comparison between theory and experiments is mostly 
qualitative. Overall, we are satisfied that the global be- 
haviour, like the increase of current with illumination, and 
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the shape of experimental curves, is weU reproduced. In 
details, figs. 3(a), 3(b) and 3(c) report the comparison of 
the I-V characteristics of bR/PC junctions [3], obtained 
by a fine tuning of the barrier height with, respectively, 
$= 53, 59, 69 meV. For $= 59 meV (Fig. 3(b)) calcula- 
tions well reproduce the /- V characteristic of the native 
state, but underestimate that of the activated state. On 
the other hand, a further increase to <i> = 69 meV does not 
fully account for the super linear behaviour (see fig. 3 (c)). 
Therefore, in analogy with the case of disordered organic 
materials [7,16], fig. 3 (d) reports the results obtained 
by taking a Gaussian distribution of $ within an average 
value of 69 meV and a dispersion cr — 44 meV. We con- 
clude that the results reported in fig. 3 (in particular fig. 
3 (d)) capture the essence of the experiments and validate 
the conjecture that the conformational change of bacte- 
riorhodopsin can be reliably detected by current transport 
measurements in the presence and absence of light. 

Remarkably, we notice that, when scaled by a factor 
of 25, the values used here for the electron effective mass 
and barrier heights are in reasonable agreement with the 
values obtained within a simple direct tunneling model 
of charge transport like in a metal-insulator-metal struc- 
ture [5] . Indeed, the factor 25 represents the mean number 
of sequential steps made by a carrier when going from one 
contact to the opposite one. Clearly, the advantage of the 
present model stems from the strict correlation between 
the microscopic structure of the protein and its macro- 
scopic electrical property. 

By construction, the present model reproduces also 
other details of experiments summarized in the introduc- 
tion. 

Theory predicts that in absence of a conformational 
change the protein resistance remains the same. This is 
verified experimentally in several ways: i) By substituting 
the original chromophore with one that is not sensible to 
the light, and finding that the I-V sensitivity to the green 
light is washed out. ii) By substituting the retinal with 
indeed with retinaloxime [3,4]. Then green light does not 
modify the conformation and thus does not produce vari- 
ation in the (nonlinear) current, iii) By using an artificial 
pigment derived from "locked" retinal. In this case the 
protein reaches an M state different from the wild M state 
and there is only a small current increase [4] . 

Theory predicts the presence or absence of the protein 
conductance on the basis of the range of interaction be- 
tween amino acids. Absence of conduction is achieved by 
taking a small interaction radius (typically less than 3.8 
A ) . This is verified by the experiments carried out on pro- 
teins deprived of the retinal, where the response to an ex- 
ternal bias is found to be a very low and completely noisy 
signal [3,4]. The origin of this result is interpreted as the 
protein loss of connectivity i. e. as a dramatic deterioration 
of its structure. 

Theory predicts the symmetry of the /- V characteris- 
tic in sweeping from negative-to-positive bias, which is in 
accordance with experimental results. 



To further test the physical plausibility of the model 
presented here, we have considered the case of (bovine) 
rhodopsin, a photoreceptor pertaining to the G-protein- 
coupled receptor (GPCR) family. Although the mecha- 
nism of photo-activation is essentially different from that 
of bR, this protein shares with bR a quite similar tertiary 
structure, i.e. both are seven-helices transmembrane pro- 
teins. This fact induced the researchers to take the bR 
structure as a template for reconstructing the 3D confor- 
mation of rhodopsin (and more generally of all the GPCR 
proteins). Accordingly, calculations have been carried out 
for the native and activated state of rhodopsin, whose 
structures were engineered as reported in [14]. The ob- 
tained results, which parallel those of fig. 3(b), are re- 
ported in fig. 4. Remarkably, for the analogous conforma- 
tional change, the model predicts that rhodopsin exhibits 
a behaviour opposite to that of bR, i. e. in going from the 
natural to the activated configuration, current at a given 
voltage is suppressed. Physically, we believe that this op- 
posite behaviour follows from the opposite change the reti- 
nal undergoes when absorbing a photon. In rhodopsin the 
retinal shape goes from bent to straight, while the reverse 
occurs in bR [12]. As a consequence, we arrive at the plau- 
sible conclusion that an opposite conformational change of 
the entire protein should imply an opposite change of its 
I-V characteristics. 

Conclusions. — We have presented a resistor network 
model implemented for a sequential tunneling mechanism 
to investigate charge transport, and in particular its struc- 
ture dependence, in proteins. The model proves its ca- 
pability to capture the correlation between a change of 
the protein conformation and that of its current-voltage 
characteristic. The theory is validated on recent exper- 
iments performed on a nanojunction filled with bacteri- 
orhodopsin, a light- activated protein, which evidences a 
substantial increase of the current, at a given voltage, in 
going from dark to the presence of a green light. A good 
agreement with experimental data is obtained when con- 
sidering an interaction radius between the amino acids of 
6 A and an average barrier energy of 69 meV with a Gaus- 
sian distribution of standard deviation cr = 44 meV. Inter- 
estingly, the model is able to explain other details of exper- 
iments, such as the role played by the presence or absence 
of the chromophore responsible of the sensing activity, as 
well as it is compatible with simple direct tunnelling mod- 
els developed on the metal-insulator-metal scheme. Under 
analogous conditions, application of the model to the case 
of bovine rhodopsin predicts an opposite behaviour with 
a suppression of the current. Confirmation of such a pre- 
diction remains an experimental challenge. 
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Table 1: Bacteriorhodopsin PDB entries 



Entry 


State 


Resolution 


(A) 


Temperature (K) 


Amino acids 


IFBB 
IFBK 


Native 
Activated 


3.20 
3.20 




— 


224 
225 


1QM8 
IDZE 


Native 
M State 


2.50 
2.50 




100 
100 


229 
224 


2NTU 
2NTW 


Native 
L State 


1.53 
1.53 




100 
100 


222 
222 


IMOK 
IMOM 


K State 
M State 


1.43 
1.43 




100 
100 


222 
224 
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Fig. 1: Running average of the network resistance for $ = 59 
meV. Data are obtained with an apphed voltage of 0.1 and 1 
V, respectively. 



factory Neuron Device (BOND) project. 
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Fig. 3: I-V characteristics of the native and activated state 
of bR, 2NTU and 2NTW, respectively. Full squares refer to 
the native state and open circles to activated state; dashed 
curves are guides to the eyes, full curves refer to experiments 
[3]. Calculations refer to different barrier heights: (a) $ = 53 
meV, (b) $ = 59 meV, (c) $ = 69 meV, (d) average $ = 69 
meV Gaussian distributed with a variance cr = 44 meV. All 
the currents obtained by calculations are normalized to the 
experimental value of current in the dark at 1 V [3]. 
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Fig. 4: Predicted /- V characteristics of the native (full squares) 
and activated (open circles) state of bovine rhodopsin, dashed 
curves are guides to the eyes. Calculations are carried out 
analogously to those in fig. 3(b). All the currents obtained 
by calculations are normalized to the experimental value of 
current in the dark at 1 V [3]. 
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